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j«-^ ■ Abstract Ground Level Enhancement (GLE) events represent the most energetic class 

^sj ' of solar energetic particle (SEP) events, requiring acceleration processes to boost ^ 1 

GeV ions in order to produce showers of secondary particles in the Earth's atmosphere 

J with sufficient intensity to be detected by ground-level neutron monitors, above the 

Mm , background of cosmics rays. Although the association of GLE events with both solar 

^^ ' flares and coronal mass ejections (CMEs) is undisputed, the question arises about the 

f~| , location of the responsible acceleration site: coronal flare reconnection sites, coronal 

CL(' CME shocks, or interplanetary shocks? To investigate the first possibility we explore 

JL I the timing of GLE events with respect to hard X-ray production in solar flares, con- 

5_( ■ sidering the height and magnetic topology of flares, the role of extended acceleration, 

"v^ ' and particle trapping. We find that 50% (6 out of 12) of recent (non-occulted) GLE 

C^ , events are accelerated during the impulsive flare phase, while the remaining half are 

accelerated significantly later. It appears that the prompt GLE component, which is 

j^-v , observed in virtually all GLE events according to a recent study by Vashenyuk et al., 

^ ' is consistent with a flare origin in the lower corona, while the delayed gradual GLE 

^\ ' component can be produced by both, either by extended acceleration and/or trapping 

Cn ■ in flare sites, or by particles accelerated in coronal and interplanetary shocks. 

O 

c7^ ' Keywords Solar Flares — Particle Acceleration — Ground Level Enhancements 

in 

o. 

^— ^ ■ 1 Introduction 



A key aspect that motivated this review is the question whether ground level enhance- 
ment (GLE) events, which apparently require acceleration processes that produce ^ 1 
, ■ , GeV particles, originate from flare regions in the solar corona or from shocks driven by 

j^ ' coronal mass ejections propagating through the corona and interplanetary space. GLE 
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events represent the largest solar energetic particle (SEP) events that accelerate GeV 
ions with sufficient intensity so that secondary particles are detected by ground-level 
neutron monitors above the galactic cosmic-ray background (Lopate 2006; Reames 
2009b). A catalog of 70 GLE events, occurring during the last six solar cycles from 
1942 to 2006, has been compiled (Oliver et al. 1982; Oliver 2006), which serves as the 
primary database of many GLE studies. So, GLE events are very rare, occurring only 
about a dozen times per solar cycle, which averages to about one event per year. While 
GLE events with 1 GeV energies represent the largest energies produced inside our so- 
lar system, they are at the bottom of the cosmic ray spectrum, which covers an energy 
range of « 10 — 10 eV, exhibiting a "spectral knee" between particles accelerated 
inside our galaxy (« 10 — 10 eV) and in extragalactic sources (~ 10 — 10 eV). 
While coronal mass ejections (CMEs) are widely considered as the main drivers of 
geoeffective phenomena, as pointed out in the so-called "solar flare myth" paradigm 
(Gosling 1993), the acceleration site of high-energy particles detected in-situ in the he- 
liosphere can often not unambiguously be localized, and thus we have to consider both 
options. Several timing studies of GLE events have shown evidence for multiple (impul- 
sive and gradual) components, so there is often not a simple dichotomy of acceleration 
sites. In order to address our central question whether GeV particles producing GLE 
events are either of solar or of heliospheric origin, we will review various aspects of flare 
observations of GLE events, such as the relative timing of the release of OLE-associated 
particles during solar flares (Section 2.1), prompt flare-associated acceleration of GLE 
protons (Section 2.2), the maximum energies of gamma-ray producing particles in so- 
lar spectra (Section 2.3), the height of acceleration regions (Section 2.4), the magnetic 
topology of solar flare acceleration regions (Section 2.5), and the role of extended ac- 
celeration phases and particle trapping in solar flares (Section 2.6). We summarize 
our conclusions in Section 3. Complementary aspects of the same question from the 
view of CME-associated shocks are treated in the article of Gang Li, and active region 
characteristics of GLE-associated flares are reviewed by Nariaki Nitta in this volume. 



2 Flare Observations of GLE Events 

All GLE events are associated with solar flares of the most intense category, i.e., GOES 
X-class flares in most cases, althouth there are exceptions, e.g., see the 1979 August 
21 event (Oliver et al. 1983) or the 1981 May 10 event (Oliver 2006). At the same 
time, coronal mass ejections (OME) were reported in all recent cases. Thus we can say 
that flares and OMEs are both necessary conditions for a GLE event, but it leaves us 
with the ambiguity where the acceleration of GeV particles reponsible for GLE events 
takes place. In the following we investigate and review various observational aspects of 
relevant flare data that could shed some light into this question. 



2.1 Solar Particle Release Times 

For the analysis of the timing of GLE events we are using initial solar release times 
that are normalized to the Earth-arrival times of electromagnetic signals, based on the 
compilation of 16 events during 1997-2006 (Gopalswamy et al. 2010; Reames 2009a) 
and of 30 events during 1973-2006 (Reames 2009b). 
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Fig. 1 Time profiles of tfie GOES 1-8 A soft X-ray fluxes (solid red curve), GOES 0.5-4 A 
fluxes (dashed red curve), and their time derivatives (thick black curves), which are a proxi 
for the hard X-ray flux, and solar particle release times with uncertainties (hatched areas) 
of 13 GLE events during 1998-2006. The time range is centered at the hard X-ray peak time 
(given on the left side) . The soiar particle release times (SPR) are determined from the velocity 
dispersion of particles detected with WIND by Reames (2009b), corrected by -1-8.3 min to the 
photon arrival time at Earth. The start times of radio type II bursts are indicated with vertical 
dashed bars. 



Table 1 Timing of 13 GLE events during 1998-2006 (corrected to photon arrival time at 
Earth). 



GLE 


Date 


Type II 


Hard X-ray 


SPR 


Delay 


Delay 


Location 


# 




start time 


peak time 


time 


SPR-II 


SPR-HXR 








tii [UT] 


tHXR [UT] 


tsPB. [UT] 


[min] 


[min] 




56 


1998-May-02 


13:41 


13:39:39 


13:55:02 


4-14.0 


+15.4 


S15W15 


57 


1998-May-06 


08;03 


08:05:12 


08:11:50 


+8.8 


+6.6 


S15W64 


58 


1998-Aug-24 


22:02 


22:05:54 


22:40:26 


-1-38.4 


+34.5 


N35E09 


59 


2000-Jul-14 


10:19 (10:28*) 


10:18:45 


10:24:50 


-3.2 


+6.0 


N22W07 


60 


2001-Apr-15 


13:47 


13:47:06 


13:54:02 


-H7.0 


+8.9 


S20W84 


61 


2001-Apr-18 


02:17 


02:13:27 


02:32:38 


-1-15.6 


+19.2 


«W115 


63 


2001-DCC-26 


05:12 


05:17:30 


05:28:56 


-1-16.9 


+11.4 


N08W54 


64 


2002-Aug-24 


01:01 


01:02:33 


01:08:26 


+7.4 


+5.8 


S02W81 


65 


2003-Oct-28 


11:02 


11:05:12 


11:13:26 


+11.4 


+8.2 


S16E08 


66 


2003-Oct-29 


20:42 


20:44:27 


21:03:56 


+21.9 


+19.5 


S15W02 


67 


2003-NOV-02 


17:14 


17:19:00 


17:22:08 


+8.1 


+3.1 


S14W56 


69 


2005-Jan-20 


06:44 


06:48:27 


06:47:50 


+3.8 


-0.7 


N12W58 


70 


2006-Dec-13 


02:26 


02:27:24 


02:42:20 


+16.3 


+14.9 


S06W26 



References: Reames et al. 2001 (#59); Aschwanden and Alexander 2001 (#59); Pohjolainen 
et al. 2001 (#56); Shumilov et al. 2003 (#56); Vashenyk et al. 2003 (#59); Chertok et 
al. 2001 (#59); Torsti et al. 2004 (#56); Miroshnichenko 2005 (#66, 67); Klassen et al. 2005 
(#65); Simnett an Roelof 2005a, b (#69); Kocharov et al. 2006 (#56); Martirosyan and 
Chilingarian 2005 (#69); Kuznetsov et al. 2006 (#69); Simnett 2006, 2007 (#69); Wang 
and Wang 2006 (#69); Li et al. 2007 (#59); Bieber et al. 2008 (#69,70); Bombardieri et 
al. 2008 (#69); Grechnev et al. 2008 (#69); McCracken et al. 2008 (#69); Moraal et al. 
2008 (#69); Watanabe et al. 2008 (#65); Chupp and Ryan 2009 (#65,69); Li et al. 2009 
(#70); Masson and Klein 2009 (#69); Matthid et al. 2009 (#69); Reames et al. 2009a,b (all 
events); Wang 2009 (#59); Vashenyuk et al. 2011 (#59-61,65,67,69,70). - *) There are two 
type II events for the 2000-Jul-14 flare, one starting at 10:19 UT (Chertok et al. 2001), the 
other starting at 10:28 UT (Gopalswamy 2010). 



In Fig. 1 we display the time profiles of 13 events (out of the 16 GLE events 
occurring during 1994-2007) analyzed in Reames (2009a), for which the solar release 
time could be determined with the velocity dispersion method. The flare time profiles 
we show in Fig. 1 are the soft X-ray fluxes IsxRif) detected by GOES in the soft (1-8 
A) and hard (0.5-4 A) channel (red solid and dashed curves), as well as their time 
derivatives (black solid and dashed curves), which are generally considered as a good 
proxi of the hard X-ray flux InXRit), 



iHXRit) 



dlsXRJt) 
dt 



(1) 



according to the Neupert effect (Neupert 1968; Hudson 1991; Dennis and Zarro 1993; 
see also Section 13.5.5 in Aschwanden 2004 and references therein). Note that the time 
derivatives of both soft X-ray channels yield a near-simultaneous hard X-ray peak time 
(< 1 min) within the accuracy needed here. The time profiles shown in Fig. 1 are all 
of the same length (1 hr) and centered at the hard X-ray peak time tfjxRy which we 
consider as the reference time for acceleration of the most energetic electrons during 
a flare. The physical explanation for the Neupert effect is that the hard X-ray peak 
time coincides with the time of most intense precipitation of accelerated high-energy 
electrons down to the chromosphere (within a time-of-flight interval of « 10— 100 ms), 
which represents the time interval of most intense heating of chromospheric plasma, 
giving rise to the steepest increase of the soft X-ray flux during a flare. These hard 
X-ray peak times t^xR a-re hsted in Table 1, which typically occur about ~ 5 — 10 
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Fig. 2 Example of a velocity dispersion measurement of GLE arrival times of protons at 
IMP-7, IMP-8, or neutron monitors (y-axis in left panel) versus the reciprocal velocity (l/v) 
(x-axis in left plot), which allows to extrapolate the energy release time tspR at the source 
location (1/v i— >■ 0), for the GLE event of 1978 May 7. The timing of the relative intensity of 
the protons with different energies is shown in the right-hand-side panel (Reames 2009b). 



min after the soft X-ray onset times, but preceed the soft X-ray peak times by about 
1-15 min. The start or peak of soft X-ray emission (as listed in GOES flare catalogs) 
should not be used as a reference time for particle acceleration, because they rather 
bracket the beginning and end of significant hard X-ray emission and the concomitant 
acceleration phase. 

The solar particle release times tgpf; on the other hand, were measured from 
the velocity dispersion of particles arriving at ground-based neutron monitors or at 
particle detectors in near-Earth satellites, such as with the Interplanetary Monitoring 
Platform 8 (IMPS) and the WIND spacecraft (Reames 2009a, b) . An example of such a 
measurement is shown in Fig. 2 for the GLE event of 1978 May 7. A plot of the arrival 
time versus the reciprocal velocity (l/v) shows a nearly-linear relationship that can be 
fitted with a linear regression fit and yields the path length (~ 1.1 — 2.1 AU) as well as 
the so-called solar particle release time tgpn within an accuracy of typically ±1, ..., 3 
min. We list these extrapolated solar particle release times tspR determined by Reames 
(2009a) in Table 1 and mark their range of uncertainty with a hatched area in Fig. 1. 
Defining the time interval of the impulsive flare phase where the time derivative of 
the soft X-ray flux is positive (corresponding to the time interval of impulsive hard 
X-ray emission according to the Neupert effect), we flnd that 6 out of the 12 events 
(50%) have the earliest particle release times tspp overlapping with the impulsive 
hard X-ray flare phase (ignoring the occulted event #61 on 2001-Apr-18), and thus are 
consistent with flare-associated acceleration, which is generally called the prompt GLE 
component. The remaining events have particle release times that start between 8 and 
34 minutes later than the hard X-ray peak, which could be produced by acceleration 
in CMEs, based on the relative timing to hard X-rays. This results is not inconsistent 
with an earlier review by Oliver et al. (1982), where the most likely injection onset of 
GeV protons was found to conincide with the flrst significant microwave peak (which 



is produced by gyrosynchrotron emission and generally coincides with the production 
of nonthermal hard X-rays) . 

Another timing indicator is the onset time tu of metric type II bursts, which is 
also listed in Table I (from Gopalswamy et al. 2010). The onset of type II bursts 
indicate the formation of a CME-related shock front. These type II onset times t// are 
indicated with a dashed vertical bar in Fig. 1, which seem to occur mostly during the 
impulsive flare phase or rise time of soft X-ray emission. Apparently, shocks in these 
intense X-class flares originate promptly during the flare, about at the same time when 
the magnetic field becomes stretched out below the erupting filament, which triggers 
magnetic reconnection, particle acceleration, and chromospheric hard X-ray and soft 
X-ray emission subsequently. In half of the events, the solar particle release times tgpji 
coincide with both the impulsive hard X-ray flare phase as well as with the onset time 
tjj of type II bursts (Fig. 1), and thus it is ambiguous whether the prompt component 
of GLE particles is accelerated in the flare region or in a near-simultaneous type II or 
CME-associated shock front. 

One GLE event, for which the timing was most extensively determined is the 
second-latest event of 2005 Jan 20. The steepest increase of the soft X-ray flux oc- 
curs around 06:48 UT (Fig. 1). Hard X-ray time profiles from RHESSI show a peak of 
the 50-100 keV and 300-800 keV emission around 06:45-06:46 UT, while the gamma-ray 
emission at 2-6 MeV, 23-40 MeV, and > 60 MeV detected with SONG peak around 
06:46-06:47 UT (Kuznetsov et al. 2006; Grechnev et al. 2008). The solar particle re- 
lease time was determined at 06:39 UT at the Sun (Reames 2009a), corresponding to 
a photon arrival time of 06:47 UT at Earth, which coincides with the hard X-ray and 
gamma-ray peaks, and about 3 min after the type II start time at 06:44 UT. Based 
on the temporal and spectral properties it was concluded that the acceleration site 
leading to the SEP/GLE spike is likely to be located in the flare region rather than in 
CME shocks, at least for the leading SEP/GLE spike (Grechnev et al. 2008; Simnett 
and Roelof 2005a, 2005b; Simnett 2006, 2007; Kuznetsov et al. 2006; Wang and Wang 
2006; Bombardieri et al. 2008; Chupp and Ryan 2009; Masson and Klein 2009), while 
a second component later on could be accelerated in a CME-associated shock at a 
distance of « 3 — 5 solar radii (McCracken et al. 2008; McCracken and Moraal 2008; 
Moraal et al. 2008). 



2.2 Prompt Flare- Associated Acceleration of GLE Protons 

Most GLE events exhibit a prompt component (PC) and a delayed component (DC), 
which were identified in nearly all events in a recent study of 35 large GLE events during 
the period of 1956-2006 (Vashenyuk et al. 2011). The prompt component prevails at the 
beginning of the event and is characterized by an impulsive profile, strong anisotropy, 
and by an exponential energy spectrum, i.e. J{E) oc exp(— _B/i5o) with Eq ~ 0.5 GeV 
(within a range of 0.3 GeV < Eq <1.8 GeV). The delayed component dominates during 
the maximum and decay phase of the events, has a gradual intensity profile, a moderate 
anisotropy, and a powerlaw energy spectrum (with a typical slope of 5 « 5 ± 1). In 
comparison, the study of Moraal and McCracken (2011) identifies a first (prompt) 
component in 14 events and a secondary (delayed) component in 6 (out of 16 GLE) 
events during Cycle 23. There is discrepancy in 3 events, where Vashenyuk et al. (2011) 
registered a delayed component, while Moraal and McCracken (2011) do not detect it. 



Since CME-associated shocks last much later than the impulsive flare phase, shock- 
accelerated particles are likely to increase in number and are subject to a gradual release 
as long as the shock lasts, and thus are likely not to have an impulsive time profile, 
while flare-associated hard X-rays exhibit the same impulsive time profile of particle 
acceleration naturally. However, some theoretical models of CME-driven shocks that 
include the magnetic field gradient and the evolution of the shock-normal angle from 
quasi-perpendicular to quasi-parallel orientation predict impulsive time profiles also 
(e.g., Zank, Rice, and Wu 2000). Nevertheless, the fact that most GLE events (29 out 
of 35) analyzed in Vashenyuk et al. (2011) exhibit a prompt component, togehter with 
our finding that the GLE start times occur during the impulsive hard X-ray phase 
in 50%, is consistent with the interpretation of flare-associated acceleration for the 
prompt component. 



2.3 Maximum Energies of Solar Gamma Rays 

There is no known high-energy cutoff of the electron bremsstrahlung spectrum; the 
highest energies of observed bremsstrahlung are around several 100 MeV (Forrest et 
al. 1985; Akimov et al. 1991, 1994a,b,c, 1996; Reames et al. 1992; Dingus et al. 1994; 
Trottet 1994; Kurt et al. 1996; Rank et al. 2001), see also reviews by Ramaty & 
Mandzhavidze (1994) or Chupp and Ryan (2009). Gamma-rays were reported up to 
energies above 1 GeV with the Energetic Gamma-ray Experiment Telescope (EGRET) 
on CGRO during the GLE event #51 on 1991 Jun 11 (Kanbach et al. 1992). The 
spectrum of the flare (Fig. 3) could be fitted with a composite of a proton generated 
pion neutral spectrum and a primary electron bremsstrahlung component (Kanbach et 
al. 1992). The GLE event #69 on 2005 Jan 20 reported gamma rays, protons, and pion 
decay radiation up to energies > 200 MeV (Grechnev et al. 2008). In Table 2 we list the 
energy ranges of observed photons in hard X-rays and gamma rays and the required 
primary particle energies for the various high-energy processes operating in large solar 
flares (adapted from Ramaty and Mandzhavidze 1994). While we have observational 
information of maximum energies of electrons above 1 GeV, the maximum energies of 
protons and ions in this regime can only be estimated indirectly, e.g., by the energy equi- 
partition argument, or by numerical simulations of a theoretical acceleration model. 

The prompt component of GLE proton spectra were found to extend to maxi- 
mum energies in the range of 0.3 GeV < Eq < 1.8 GeV (Vashenyuk et al. 2011). 
Ground-based neutron detectors detect GLE protons at energies E ^ 1.5 GeV, which 
corresponds to a Lorentz factor of 7 = E/{mpC ) = 1.6, a kinetic energy of i5j.j„ = 
mpc^{^ - 1) = 562 MeV, and relativistic speed of /? = ^(1 - 1/7^) ^ 0.78. 

Do these highest observed energies constrain or rule out any acceleration mecha- 
nism? For DC electric fleld acceleration in sub-Dreicer fields, the maximum velocity to 
which electrons can be accelerated is limited by the value of the Dreicer field, which 
depends on the density and temperature of the plasma, Ejj ~ 2 x 10~ n/T (stat- 
volts cm~ ). Holman (1996) argues that electron energies up to 10 — 100 MeV can 
be attained for high densities of He ~ 10 cm~ and low temperatures T ~ 2 MK 
(yielding a Dreicer field of Ejj = 1 x 10~ statvolt cm~ , i.e., 3 V m~ ), if electrons 
are continuously accelerated over a current channel with a length of i = 10— 100 Mm. 
The requirement for such large-scale DC electric fields, however, conflicts with the ob- 
served time-of-flight delays of hard X-ray pulses (e.g., Aschwanden et al. 1996) and 
the short inductive switch on/off time scales required for the observed subsecond hard 



Table 2 Hard X-ray and gamma-ray emission mechanisms operating in large solar flares, 
with observed photon energy and primary particle energy ranges (adapted from Ramaty and 
Mandzhavidze 1994) 



Process 



Observed photon 



Primary particle 
energies 



Bremsstrahlung continuum 20 keV-11 MeV 20 keV-1 MeV 

>10 MeV 10 MeV-1 GeV 

Nuclear de-excitation lines 0.4,..., 6.1 MeV 1-100 MeV/nucl. 

Neutron capture line 2.2 MeV 1-100 MeV/nucl. 

Positron annihilation radiation 0.511 MeV 1-100 MeV/nucl. 

Pion decay radiation 10 MeV-3 GeV 0.2-5 GeV 

Neutrons induced in atmospheric cascades 0.1-10 GeV 0.1-10 GeV 

Neutron decay protons in space 20-200 McV 20-400 MeV 
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Fig. 3 A gamma-ray spectrum observed wfith EGRET/CGRO during the GLE event #51 
on 1991 Jun 11, 02:04 UT flare, accumulated during 03:26-06:00 UT (Kanbach et al. 1993). 
The spectrum is fitted with a combination of primary electron bremsstrahlung and pion-decay 
radiation. Note that pion decay is dominant at energies ^ 40 McV (Mandzhavidze & Ramaty 
1992). 



X-ray pulses. Alternatively, Litvinenko (1996) envisions super-Dreicer electric fields, 
which can generate arbitrarily high maximum electron energies within much smaller 
spatial scales, and thus can be consistent with the time-of-flight delays of the observed 
hard X-ray pulses. The maximum energy for relativistic electrons obtained over an 
acceleration time t is approximately. 



e(t) = eEl{t) ^ eEct 



(5) 



Litvinenko (2003) identifies plausible physical conditions with super-Dreicer fields of 
order E ^ 100 V m~ in reconnecting current sheets that lead to electron acceleration 
with gamma-ray energies of a few 10 MeV in electron-rich flares or to the generation of 
protons with energies up to several GeV in large gradual flares. A similar acceleration 
model was applied to the GLE event #59 on 2000 Jul 14 (Li et al. 2007). Also stochastic 
acceleration can generate 10 MeV electrons and 1 GeV protons (Miller et al. 1997), if 



a sufficiently high wave turbulence level is assumed (which, however, cannot easily be 
constrained by observations). Thus, the maximum observed gamma- ray energies imply 
only constraints for the acceleration mechanism of sub-Dreicer electric DC fields, but 
not for super-Dreicer DC electric field, stochastic, or shock acceleration mechanisms. 



2.4 Height of Acceleration Region 

Since we have a temporal coincidence of CLE particle acceleration with respect to 
flare hard X-ray emission in at least 50%, we turn now to the question of the spatial 
localization of acceleration sources. 

One method to estimate the height of the acceleration region of CLE particles 
(Reames 2009a,b) is based on an assumed height hji of the start of radio type II 
bursts, corrected for the propagation delay of the CME shock front with speed vcme 
during the time interval At = tgpfi — tjj between the start tjj of the radio type II 
burst and the solar particle release time tgpji, 

hsPR = hjj + vcme{^SPR — ^II) . (2) 

which is always higher than hjj, because vcme > and tgpji > t//. Reames (2009a,b) 
assumed a nominal average height of hgpji = 1.5 ± 0.5 solar radii based on standard 
coronal density models ne{h) and statistical start frequencies of type II bursts around 
ly ~ 100 MHz (Kundu 1965), which corresponds to an electron density of rig ~ 10 
cm~ in the case of fundamental plasma emission, i.e., fpe ~ 8980y'n^. Taking the 
propagation delay of {tspR — t//) « 5 — 15 min and the CME speeds of vqme ~ 
1000 — 3000 km s~ into account, Reames (2009a,b) arrived at estimated heights of 
hspR « 2 — 5 solar radii, with a possible dependence on the heliographic longitude. 
This method implicitly assumes that particles are accelerated in CME shocks, whose 
location is entirely tied to the height and propagation speed of the type II and CME 
shock front. 

Alternatively, since half of the CLE events exhibit a starting time tgpp during the 
impulsive phase of hard X-ray emission, we can estimate the height of their acceleration 
region from the hard and soft X-ray data. The mildly relativistic electrons accelerated 
in a flare exhibit a time-of-flight delay between their coronal acceleration site and the 
chromospheric target region where bright hard X-ray emission is observed, depending 
on their kinetic energy. From the velocity dispersion of these energy-dependent hard 
X-ray time delays, the propagation distance and height of the acceleration region can 
be calculated, with proper correction for the geometry of the trajectory and for the 
pitch angles of the particles. 

For a CLE precursor flare, which occurred on 1992 June 25, 17:32 UT (hard X- 
ray start time) at heliographic position N10/W75, such a time-of-flight measurement 
is available (Fig. 4), based on Yohkoh/HXT and SXT observations (Aschwanden et 
al. 1996). This flare is of the GOES M1.4-class and peaked in soft X-rays at 17:54 UT. 
Concomitant hard X-ray and radio emission is also studied in Wang et al. (1995), find- 
ing a spatial separation of 35 Mm between the simultaneous hard X-ray and microwave 
emission, which provides also an approximate scale for the horizontal and vertical ex- 
tent of the flare region. This flare is a precursor to a CLE event 2 hrs later at the 
same location, which is listed as CLE event # 53 in Oliver (2006), peaking at 20:14 
UT at position N10/W68, and is classified as X3 GOES-class flare. We can consider 
the two rapidly following events as homologous Bares within the same active region 
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Fig. 4 Altitude measurement of the acceleration source during the 1992 June 25 GLE precur- 
sor flare based on the velocity dispersion of electron timc-of-fiight delays, measured with the 
Compton Gamma Ray Observatory (CGRO), Yohkoh/HXT, and SXT. A hard X-ray image 
from HXT/Ml (23-33 keV) of the flare footpoint sources (contours) and a SXT soft X-ray im- 
age from Yohkoh/SXT (greyscale) are shown in the top left panel. The energy-dependent hard 
X-ray time delays in the range of At ^ 300 ms are shown in the bottom left panel, yielding 
a time-of-flight distance of L' = 58.5 ± 14.0 Mm, which is projected onto the inferred loop 
geometry (top right panel). The uncertainty of the TOF distance is derived from varying the 
fllter time scale (bottom right panel), (Aschwanden ct al. 1996). 



that are likely to have a similar magnetic configuration. The time-of-flight distance 
for the precursor flare was evaluated from the time delay r^ between electrons with 



relativistic speeds of , 
keV, 



and /3j in the hard X-ray photon energy range of e ~ 20 — 80 



TOF 



1 

To 



(3) 



which after correction for electron pitch angle, helical twist of field lines, and projection 
effects of loop size yielded a value of L' = 58.5 ± 14.0 Mm, corresponding to a height 
of ft ~ 50 ± 13 Mm. This is a fairly typical height of the acceleration region for large 
flares, amounting to about the double height of the soft X-ray loops. From statistics 
of 42 flares, an average height ratio of h/hig^p ~ L/Lig^p = 1.4 ± 0.3 was obtained 
(Aschwanden et al. 1996), for flare loop radii of ri^gp ~ 2 — 20 Mm. Thus, the height 
range of acceleration regions in flares amounts to ft ~ 4 — 40 Mm, which corresponds to 
^ 5% of a solar radius. In summary, since about 50% of the GLE events are consistent 
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Fig. 5 Top: tracings of individual flare loops from TRACE 171 A images of the Bastille-Day 
flare 2000-Jul-14. The five sets of loops traced at five different times are marked with different 
greytones. Note the evolution from highly sheared to less sheared loops. Bottom: the position 
of the two flare ribbons traced from 171 A images. Note the increasing footpoint separation 
with time (Aschwanden 2002). 



with a particle release time during the flare hard X-ray phase, they are expected to 
have acceleration heights of ft J: 0.05 solar radii. 



2.5 Magnetic Topology of Acceleration Regions 



The geometry and magnetic topology of a solar flare region has been studied in most 
detail for the GLE event of 2000 July 14, 10:18 UT, an X5.7 GOES-class flare (e.g., 
Aulanier et al. 2000; Aschwanden and Alexander 2001; Aschwanden 2002; Yan and 
Huang 2003; Aschwanden and Aschwanden 2008). The solar energy release times 
of the particles responsible for the GLE spike during the Bastille-Day flare is well- 
synchronized with the rise phase of hard X-ray emission (Fig. 1, Table 1; Reames 
2009b; Li et al. 2007; Wang 2009). The Bastille-Day flare was observed near disk cen- 
ter, which provided an excellent view on the projected flare area, which consists of 
a classical double-ribbon structure as it is typical for most large flares. The double- 
ribbon structure straddles along a neutral line from east to west, and spreads apart as a 
function of time (Fig. 5, bottom), as expected in the Carmichael- Sturrock- Hirayama- 
Kopp- Pneuman (CSHKP) standard flare model. While the CSHKP model essentially 
describes the vertical evolution in a 2-D cross-section of the flare arcade, the Bastille- 
Day flare in addition shows also the horizontal projection and the 3-D evolution in 
detail. From loop tracings at different times during the flare (Fig. 5, top) it is evident 
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Fig. 6 Scenario of the dynamic evolution during the Bastille-Day 2000-Jul-14 flare: (a) low- 
lying, highly sheared loops above the neutral line become unstable first; (b) after loss of 
magnetic equilibrium the filament jumps upward and forms a current sheet according to the 
model by Forbes and Priest (1995). When the current sheet becomes stretched, magnetic 
islands form and coalescence of islands occurs at locations of enhanced resistivity, initiating 
particle acceleration and plasma heating; (c) the lowest lying loops relax after reconnection and 
become filled due to chromospheric evaporation (loops with thick linestyle); (d) reconnection 
proceeds upward and involves higher lying, less-sheared loops; (e) the arcade gradually fills up 
with hot fiare loops; (f) the last reconnecting loops have no shear and are oriented perpendicular 
to the neutral line. At some point the filament disconnects completely from the fiare arcade 
and escapes into interplanetary space (Aschwanden 2002). 
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that low-lying, highly-sheared loops over the neutral line brighten first, triggering a 
sequence of flare loops that progresses to higher-lying and less-sheared field lines above 
the neutral line, until we see a final double-ribbon flare arcade at orthogonal angles to 
the neutral line, spanning over a width of w ~ 50 Mm in NS direction and a length 
of I ~ 250 Mm in EW direction. From this time evolution we can reconstruct the 3-D 
geometry of the magnetic field lines that are involved in the flare as shown in Fig. 6. 
Since the bright EUV postflare loops outline the relaxed fleld line configuration after 
magnetic rcconnection, we can directly infer the size and height of the magnetic recon- 
nection region, which, statistically, is located about a factor of 1.5 above the soft X-ray 
and EUV postfiare loops (Aschwanden et al. 1996), estimated to be in an altitude of 
h ^ w ~ 50 Mm with a length of I ~ 250 Mm. This height range of the magnetic rccon- 
nection region confines the overall volume of the particle acceleration region, although 
a highly fragmentated structure with many magnetic islands is expected, caused by 
the tearing- mode instability and bursty reconnection (e.g., Sturrock 1966; Karpen et 
al. 1995, 1998; Kliem et al. 2000; Shibata and Tanuma 2001). 

Some consequences of this magnetic topology for GLE events are: (1) The vertical 
X-type reconnection geometry allows particle acceleration in upward and downward 
direction in a quasi-symmetric fashion, so that particles of almost equal energies can 
be accelerated in both directions; (2) The magnetic field lines above the vertical current 
sheet of the main reconnection regions are likely to be open, which allows escape of 
accelerated particles into interplanetary space and along Earth-connected magnetic 
field lines; (3) The relatively low height in the solar corona {h J; 50 Mm) provides a 
large reservoir of thermal particles that can be accelerated in sufficient numbers that 
are necessary for GLE detection. The escape conditions of accelerated particles into 
interplanetary space requires open magnetic field lines, which exist not only in coronal 
holes but also to a substantial fraction in active regions. Schrijver and DeRosa (2003) 
found from potential-field extrapolations of the global magnetic field over the entire 
solar surface that the interplanetary magnetic field (IMF) originates typically in a dozen 
disjoint regions, around the solar cycle maximum. While active regions are often ignored 
as a source for the interplanetary magnetic field, Schrijver and DeRosa (2003) found 
that the fraction of the IMF that connects directly to magnetic plages of active regions 
increases from ^ 10% at cycle minimum up to 30 — 50% at cycle maximum, with even 
direct connections between sunspots and the heliosphere. Evidence for open-field escape 
routes was demonstrated with magnetic field modeling, for instance for GLE event #70, 
2006 Dec 13, 02:27 UT (Li et al. 2009). In flaring regions that expel a CME, particles 
accelerated in the reconnection region behind the rising filament may be trapped inside 
the CME bubble and cannot directly escape into interplanetary space (Reames 2002), 
but may contribute to a seed population for secondary acceleration in CME-associated 
shocks (McCracken and Moraal 2008). The two populations of directly-escaping and 
trapped-plus-accelerated particles may be distinguishable as impulsive and gradual 
phases of SEP events. 



2.6 Extended Particle Acceleration and Trapping 

The solar particle release time tgpn as determined by Reames (2009b) occurs after 
the peak time tm of the hard X-ray production phase in 92% (12 out of 13 cases), or 
after the end time te of impulsive hard X-rays in 50% (6 out of 12 non-occulted cases), 
for the GLE events listed in Table 1 and shown in Fig. 1. The question arises whether 
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Fig. 7 Extended 7-ray emission as measured by COMPTEL for the GLE event #51 on 1991- 
Jun-11 fiare in the 2.223 MeV neutron capture line (top) and the 4 — 7 MeV nuclear line 
flux (bottom). The data have been corrected for primary and secondary brcmsstrahlung. A 
two-fold exponential decay has been fitted. The origin of the time axis is 01;56 UT, and the 
flare onset reported by GOES. Only data of the extended phase (after 02:13 UT) arc shown 
(Rank et al. 2001). 



we can exclude an interpretation in terms of a flare acceleration site for such delays. 
If the flare site shows, besides the impulsive hard X-ray phase, also prolonged time 
intervals with hard X-ray and gamma-ray emission, either an extended (second-step) 
acceleration phase or extended trapping is possible (e.g., Vilmer et al. 1982), in which 
case the particles responsible for GLE signatures could possibly be accelerated at the 
coronal flare site rather than in heliospheric CME shocks. Extended acceleration can be 
diagnosed from impulsive bursts in hard X-ray, gamma-ray, and radio emission after the 
impulsive phase, while indicators of extended trapping are: (1) exponential decay of X- 
ray and gamma-ray light curves, (2) spectral hardening in hard X-rays (Kiplinger 1995; 
Grayson et al. 2009), or (3) type IV continuum emission in microwaves (produced by 
gyrosynchrotron emission from trapped high-relativistic electrons). Particle trapping 
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times are limited by the collisional deflection time, which is for electrons 

UraAe) ^ t,en{e) = 0.95 ij^r^X'' f '""'"'' "! {^) (^) . (4) 



and a factor of ~ 60 longer for ions. Postflare loops usually expand slowly with height, so 
that the density drops continuously and collisional deflection times (and thus trapping 
times) become progressively longer. Particle densities in the upper corona drop below 
Tie ^ 10 cm^ , which enables trapping times over several hours. 

The second diagnostic of extended acceleration, the spectral hardening after the 
impulsive phase, has recently been demonstrated with RHESSI to be a reliable predictor 
of SEP events (Grayson et al. 2009) . From a dataset of 37 magnetically well-connected 
flares, 12 out of 18 flares with spectral soft-hard-harder (SHH) evolution produced SEP 
events, and none of the 19 flares without SHH behavior produced SEPs. Three of the 
37 analyzed SEP events are GLE events and show all SHH behavior. This demonstrates 
a statistically significant correlation between SHH and SEP (and GLE) observations. 
Since the spectral hardening happens in the hard X-rays sources at the solar flare site, 
as imaged by RHESSI, this link between SHH and SEP is unexplained in the standard 
scenario of SEP acceleration in CME-related shocks during interplanetary propagation. 

Extended gamma-ray emission with exponential decay times has indeed been ob- 
served up to 5-8 hours after the impulsive phase of the GLE event #51 on 1991 Jun 11 
(Fig. 7) and GLE event #52 on 1991 Jun 15 (Kanbach et al. 1993; Rank et al. 1996, 
2001). Evidence for prolonged acceleration of high-energy protons [E > 100 MeV) ac- 
celerated in flares and postflare loop systems, with associated spectral hardening late 
in the flare, was also discussed for GLE event #52, 1991 Jun 15, 08:17 UT (Chertok 
1995), for GLE event #55, 1997 Nov 6, 11:55 UT (Klein and Trottet 2001; Murphy 
et al. 2001; Masuda 2002; Masuda and Sato 2003), for GLE event #59, 2000 Jul 14, 
10:24 UT (Livshits and Belov 2004; Li et al. 2007; Wang 2009), and for GLE event 
#70, 2006 Dec 13, 02:27 UT (Li et al. 2009). Evidence for extended particle accel- 
eration in the lower corona was also demonstrated from microwave gyrosynchrotron 
emission, requiring high magnetic flelds near sunspots, e.g., for GLE event #69, 2005 
Jan 20, 06:48 UT (Grechnev et al. 2008; Masson et al. 2009). The GLE event #65 on 
2003 Oct 28, 11:05 UT, shows multiple phases of impulsive particle injections, ending 
with a gradual phase (> 1 hr) of type HI bursts that starts 25 min after the first 
type HI bursts (Klassen et al. 2005; Miroshnichenko et al. 2005). While there is ample 
evidence for extended acceleration of electrons in flare sites and flare-associated radio 
burst sources, the situation is less clear for ions and protons; extended acceleration of 
ions and protons can occur concomitantly in the same flare site, or alternatively in 
propagating CME shocks (e.g., Oliver 1996). 

Let us also review the four GLE events listed in Table 1 that exhibit the largest 
delays in solar particle release times tgpji. The GLE event #56 on 1998 May 02, 13:38 
UT, has a delay of tgpji — tnXR = +15.4 min, at a time when the impulsive flare 
phase is over, but the radio dynamic spectra from Artemis IV show type IV continuum 
emission at frequencies of i/ « 300 - 700 MHz during 13:42-13:50 UT (see Fig. 8 in 
Pohjolainen et al. 2001), a clear sign of gyro-synchroton radiation from trapped high- 
relativistic electrons during the GLE particle release time. 

The GLE event #58 on 1998 Aug 24, 22:05 UT, has the longest delay of tsPR ~ 
tnXR = +34.5 min, but the USAF/ RSTN data ( http://cdaw.gsfc.nasa.gov/meetings/l ws-cdaw2009/data/White/ \ 
show a late increase of 15.4 GHz radio emission after 22:30 UT, which is an indicator 
of trapped high-relativistic electrons. 
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The GLE event #61 on 2001 Apr 18, 02:12 UT, has a delay oitspR-tuxR = +19.2 
min, but this flare is occulted by 30° behind the limb, which may explain the offset 
timing due to the missed hard X-ray and soft X-ray peak. In an occulted flare we do 
not see the prompt hard X-rays from the footpoints, but only the prompt soft X-rays in 
coronal heights visible above the limb. The coronal high-temperature plasma above the 
limb cools faster (by conductive cooling) than the bulk of occulted soft X-ray plasma 
that gradually fills the postfiare loops by chromospheric evporation. Thus, the time 
profile of soft X-ray emission in occulted fiares decays faster than in unocculted flares, 
and thus we underestimate the full duration of impulsive hard X-ray emission using 
the time derivative of soft X-rays. 

The GLE event #66 on 2003 Oct 29, 20:44 UT, has a delay oitsPR-tuxR = +19.5 
min, but the USAF/RSTN data show very bursty radio spikes late in the flare during 
20:50-20:56 UT at 610 MHz, a possible indication of a secondary acceleration phase 
producing electron beams in the lower corona. In summary, we observe in all cases 
with late GLE solar particle release times signs of extended particle acceleration or 
trapping, which does not exclude that the particles responsible for the GLE signatures 
have been accelerated and temporarily trapped (for ~ 10 — 30 min) at the solar flare 
site. 



3 Conclusions 

We explored here the question whether the largest SEP and GLE events that accelerate 
ions with energies of J^ 1 GeV could be accelerated in solar flare regions, in contrast 
to the generally accepted paradigm of acceleration in heliospheric CME shocks. We 
reviewed the pro and con aspects from the solar flare site that are relevant to answer 
this question, while the complementary aspects from CME-associated shocks are dis- 
cussed in the companion article by Gang Li. The major arguments for a flare-driven 
acceleration of GeV particles are: (la) the coincidence of the particle release time with 
the impulsive flare phase (in 50%), (2a) the similarity of the impulsive time proflles 
of the flrst GLE component with the impulsive flare phase in hard X-rays, (3a) obser- 
vational evidence of maximum energies 5; 100 MeV in electron bremsstrahlung, and 
(4a) the relatively high magnetic fleld and free magnetic energy that is available in 
flare sites, compared to the locations of interplanetary CME shocks. In contrast, major 
arguments against this interpretation are: (lb) signiflcant delays of the particle release 
time with respect to the impulsive flare phase (in 50%), and (2b) the lack of direct 
observational evidence and diagnostics for ^ 1 GeV protonts and ions in flare sites. 
The conclusions are based on observations of 70 GLE events over the last six decades, 
in particular on the 13 GLE events during the last solar cycle 23 (1998-2006) that 
provided excellent new imaging data in gamma rays and hard X-rays (RHESSI), in 
soft X-rays and EUV (TRACE, SOHO/EIT), and particle data from IMP, WIND, and 
ACE. Our conclusions are: 

1. Solar particle release times tgpR of GLE-producing particles overlap with the im- 
pulsive phase of gamma-ray and hard X-ray emission in solar flares in 6 out of 12 
(non-occulted)cases, and thus the acceleration time of GLE particles is consistent 
with the flare site in 50% of the cases, taking the full duration of impulsive flare 
hard X-ray emission {tx ~ 3 — 13 min) into account. 

2. In the remaining cases, 6 our of 12 occur delayed to the flare peak by ~ 10 — 30 
min, but observational signatures of extended acceleration and/or particle trapping 
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are evident in all strongly delayed cases, and thus all GLE events could potentially 
be accelerated in flare sites. The alternative explanation of delayed second-step 
acceleration in CME-associated shocks cannot be ruled out, however, possibly con- 
stituting a secondary gradual GLE component. 

3. The height of the acceleration region of ^ 1 GeV electrons and ions depends on 
the interpretation, being h ^ 0.05 solar radii for flare site acceleration (according 
to electron time-of- flight measurements), or /i « 2 — 5 solar radii for CME shock 
acceleration (according to type II source heights and GLE solar particle release 
delays) . 

4. The magnetic topology at the particle acceleration site is not well-known from 
magnetic modeling or tracing of coronal structures, but there is statistical evidence 
for open-field regions in most active regions that provide escape routes for GLE 
particles. 

5. The recently discovered strong correlation between the spectral soft-hard-harder 
(SHH) evolution of solar hard X-rays and SEP events poses a new challenge. It is 
presently unclear how the SHH evolution can be explained in the context of the 
standard scenario in terms of SEP acceleration in CME-associated shocks. 

6. The maximum particle energies observed in solar flares reach up to several 100 
MeV for electrons and above 1 GeV for ions. The required particle energies are 

1z 1 GeV for the observed bremsstrahlung continuum at photon energies of 5; 10 
MeV, and 0.2 — 5 GeV for pion decay radiation observed at photon energies of 10 
MeV —3 GeV, which is sufficient to explain GLE detections. 

7. Energies up to < 1 GeV can be achieved for DC electric field acceleration in super- 
Dreicer fields, for stochastic acceleration, and for shock acceleration. Only DC 
electric field acceleration in sub-Dreicer fields can be ruled out, based on the large 
distances of DC fields required and the inconsistency with measured electron time- 
of-flight delays. 

What do we need to answer the question of the origin of GeV particles with more 
certainty? The velocity dispersion measurements with IMP, WIND, and ACE helped 
enormously to narrow down the start window of solar particle release times, which could 
be complemented with directivity measurements of the particle detector IMPACT on 
STEREO. The most sensitive hard X-ray and gamma-ray detectors existed on CGRO, 
but in near future we have only RHESSI available. However, RHESSI provides excel- 
lent imaging capabilities, which revealed intriguing information on displaced electron 
and ion acceleration sites in the solar corona (Hurford et al. 2003, 2006). In addition 
we expect also improved magnetic modeling of solar fiare regions that can reveal open- 
field and closed-field geometries with more certainty, using nonlinear force-free field 
(NLFFF) codes (e.g., Bobra et al. 2008; DeRosa et al. 2009) with 3-D vector magneto- 
graph data from the Solar Dynamics Observatory (SDO), possibly in conjuction with 
stereoscopic 3-D reconstruction of fiaring active regions using STEREO/EUVI. 

In conclusion, acceleration of GeV particles in fiare sites is a possibility that cannot 
be firmly ruled out with the current localization capabilities of energetic particles. Cer- 
tainly we have evidence for both acceleration in coronal fiare sites and in heliospheric 
CME shocks, often appearing concomitantly, but with different (impulsive vs. gradual) 
time scales, relative timing, and charge state characteristics. While one-sided emphasis 
has been given to both, either flares (the "big flare syndrome"), or CMEs (the "flare 
myth"; Gosling 1993), there is a consensus now that both fiare and CME phenomena 
are part of a common magnetic instability, and that both are being able to accelerate 
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particles to high energies. The remaining questions are then mostly what the relative 
proportions of both components are and how we can discriminate between them. A pre- 
liminary answer is that the observations are mostly consistent with a flare-associated 
"prompt GLE component" and a CME-associated "delayed GLE component" . 
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